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ABSTRACT 

We have studied the relation between the color evolution of Type la Supernovae (SNe la) from 
maximum light to the Lira law regime and the presence of narrow absorption features. Based on a 
nearby sample of 89 SNe la we have found that the rate of change oi B — V colors at late phases 
(between 35 and 80 days after maximum) varies significantly among different SNe la. At maximum 
light, faster Lira law B — V decliners have significantly higher equivalent widths of blended Na I Dl 
& D2 narrow absorption lines, redder colors and lower Ry reddening laws. We do not find faster 
Lira law B — V decliners to have a strong preference for younger galaxy environments, where higher 
interstellar material (ISM) column densities would be expected. We interpret these results as evidence 
for the presence of circumstellar material (CSM). The differences in colors and reddening laws found 
at maximum light are also present 55 days afterward, but unlike the colors at maximum, they show a 
significant variation among different host galaxy morphological types. This suggests that the effect of 
ISM on the colors is more apparent at late times. Finally, we discuss how the transversal expansion 
of the ejecta in an inhomogeneous ISM could mimic some of these findings. 
Subject headings: supernovae: general — distance scale 



1. INTRODUCTION 

Type la supernovae (SNe la) are thought to be the 
thermonuclear explosion of white dwarf stars. They 
are one of the main sources o f iro n in the Uni- 
verse (e.g. iMatteucci fc Greggiol I1986D and are key 
tools f or cosmologica l dista nce determinations (jPhillipsI 
Il993t iHamuv et al.1 Il996| ). leading to t he discov- 
ery of the acceleration of t he Universe (jRiess et al.l 
119981 iPerlmutter et al.l I1999 D. Despite much obser- 
vational and theoretical effort, we lack a clear un- 
derstanding about their progenitors (for review s, see 
IHillebrandt fc Niemevedr2000l: IWang fc Hanll2012D . 



Most favored progenitor scenarios involve a carbon- 
oxygen white dwarf (CO WD) in a binary system. The 
CO WD could either ignite at the Chandrasekhar mass 
after accreting matter stably from a companion in a de- 
layed detonation or pure deflagration explosion (Chan- 
drasekhar mass - single degenerate scenario, Alch - 
SD: Hachisu ct al. 1996, 1999; Hachisu, Kato. fc Nomot^ 
1999t ILi fc van den Hcuvel 1997; Langer et al. 200^ 
Han fc Podsiadlows ki 2004; Meng C hen, fc Han 200c 



stellar material left from the accretion process or from 
previous nova explosions, 2) the presence or absence of 
a companion star, 3) the SN rate delay-time distribu- 
tion, which is generally found to be inversely propor- 
tional to the time sinc e formation (for a review, see 
iMaoz fc Mannuc"cill2012l ) and is g enerally believed to fa- 
vor D D scenarios (although see iKobavashi fc Nomotol 
|2009( ). and 4) the degree of neutronization in the cen- 
tral regions due to differ ent central densities at ignition 
(e.g. 'Gera rdv et aI1l2007D . 

Although the lack of evidence for shocked compan- 
ion stars of the progenitor CO WDs, shocked presuper- 
nova winds or companion pre-explosion images is more 
generally consistent with DD scenarios (|Li et al.l 2011 



Wang et all l2009t Khokhlo^ ll99ir;Ropke & Nicmcyer 



20071 : iKromer et al.li2012i) . at a mass below the Chan- 
drasekhar mass if accretion is unstable and a surface 
detonation is strong enough to trigger a central detona- 
tion (sub-Chandra sekhar mass - s ingle degenerate sce- 
nario, Mch - SD: iNomotol 119821: iWooslev etHI 119861: 
iSim et aI1l2010t IKromer et al.ll2010l ) , or after the merger 
with a degenerate companion in a double degenerate 
scenario if the detonation can occur in a dynamical 
time -scale (violent double degenerate merger, violent 
DDHIben fc Tu tukovlgsl lWebbinlfigsllPakmor et al.l 
|2012[ ). Some observational tests for these scenarios can 
be: 1) the presence or absence of low velocity circum- 



Nugent et al."2011'; 'Bloom ct al."2012'; 'Brown et al.l l2012 ; 
Schaefer fc Pagnotta 2012; Edwards et al. 2012), there 
are different results that suggest the presence of circum- 
stellar material (CSM) around SN la progenitors, which 
is gener ally expected in S D scenarios (or slow DD merg- 
ers, see lFrver et al.|[2010r ). These consist of optical col- 
ors that favor scattering in nearby materia l (jWa ng 200 J; 
IGoobail 120081: lAmanullah fc Goobail 1201 ID. narrow ab- 
sorption lines that show variabilitv _(|Patat et al.l 120071 : 

fondin ct al. 2009; Simon ct al. 2001) are blueshifted 
ternbcrg et al.,,2011, : .Foley et al.,, 20121) . or that corre- 
late w ith certain SN ejecta viewing angles (jForster et al.l 
l2012f ). This evidence is contrasted by the absence of X- 
ray or radi o emission from sho cked CSM in the nearby 
SN 2011fe (jHoresh et al.ll2012f l. but is supported by re - 
cent observations of SN la PTFllkx (jDildav et al.ll2012| ). 
which was probably surrounded by narrow nova shells 
ejected before explosion, and the sampl e of SNe la with 
CSM interaction recently identified by [Silverman et al.l 

ipoTl . 
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If some SNe la are surrounded by significant amounts 
of CSM, it would be expected that the hght from those 
SNe would be scattered away of the line of sight by in- 
tervening dust, causing extinction and reddening, but 
also that it would be scattered into the line of sight by 
nearby dust , i.e. multiple s cattering processes or ligh t 
echoes (e.g. IChevaliedfloM IWa^[200l IGoobadl2Q08l) . 
The typical delay of these light echoes would depend on 
the typical distance to the closest surviving dust after ex- 
plosion, whi ch for ty pical distances of about 0.005 — 0.1 
pc (see e.g. iPatat e t al. 2007; Simon ct al. 2009) would 
correspond to approximately 6 to 112 days of light travel- 
time. Assuming that the light echoes and/or dust sub- 
limation at later times change the SN light curves at 
these time-scales after maximurcQ, we have done a sys- 
tematic study of the evolution o f late time colors of SN 
la du ring the Lira law regime (Lira 1995i iPhillips et al.l 
Il999f l. between 35 and 80 days after maximum. Given 
that SN colors evolve in an approximately homogeneous 
fashion during this phase, we try to detect the presence 
of multiple scattering processes, or light echoes, studying 
the evolution of colors in a big sample of nearby SNe la 
and its relation with the presence of narrow absorption 
features. 

In order to distinguish between interstellar material 
(ISM) and circumstellar material (CSM) we have stud- 
ied the relation between the equivalent widths of blended 
Na I Dl & D2 narrow absorption features, obtained from 
mid-resolution spe ctra using the method described in 
iForster et al.l 1201^ hereafter FG12; with the Lira law 
B — V decline rates. Extinction by ISM should in prin- 
ciple only shift the B — V colors to higher values, with- 
out significant changes of the Lira law B ~ V decline 
rate, whereas CSM can both shift the B ~ V colors to 
higher values and change the Lira law B — V decline rate, 
making i t generally faster as the dist ance to the CSM is 
reduced (jAmanullah fc GoobaH[20Tl or if CSM dust is 
destroyed at late times. Given that systems with signifi- 
cant CSM should have a different extinction law (jWand 
l2005t iGoo bar 2008) , we also study the relation between 
B — V and V — i colors at maximum and 55 days after 
maximum for different Lira law B — V decline rates in our 
sample. To discard that SN light curves have different 
intrinsic decline rates for different galaxy populations we 
have also performed several environmental tests studying 
the distributions of host galaxy morphological type, host 
galaxy inclination, stretch or Amis of different Lira law 
B — V decline rates. 

Energy deposition during the Lira law regime is domi- 
nated by the /?-decay of ^^Co into ^^Fe, via electron cap- 
tures (81%) or positron emissions (19%). At 35 days after 
maximum, energy deposition is dominated by the scat- 
tering of 7-rays produced in the primary channel. At 80 
days after maximum, the ejecta becomes transparent to 
7-rays and the kinetic energy from positrons produced in 
the secondary cha nnel dominates the energy deposition 
(jMilne et al.lll999() . This means that a transition from 
7-ray to positron energy deposition occurs during the 
Lira law regime, which may be linked to different bolo- 
metric decline rates. However, to our knowledge there 
is not a clear connection between the bolometric decline 

^ We define maximum as the maximum of the SN light curve in 
S— band. 



rate with the B — V decline rates observed during the 
Lira law. 

Finally, given the connection found by FG12 between 
narrow absorption features and the n ebular velocity 
shifts introduced bv iMaeda et al.l (|2010( ). interpreted as 
evidence for CSM in positive nebular velocity shift SNe 
la and a general asymmetry of the CSM post-explosion, 
we have tested whether there is any correlation between 
nebular velocity shifts and the observed Lira law B — V 
decline rates. 

2. SAMPLE SELECTION AND ANALYSIS 

We have selected those SNe la that have spectra in 
the Center for Astrophy sics (CfA) public sample of SNe 
la (jBlondin et al.l I2012D or in the Carnegie Supernova 
Prog ram (CSP) sample of nearby SNe la (Folatelli et al.l 
12013 ). For every spectrum of every SN we have run an 
automatic Python script to measure equivalent widths 
and associated errors of blended Na I Dl & D2 narrow 
absorption features as described in FG12. These values 
were combined using a weighted average to obtain one 
equivalent width and associated error per SN, assigning 
an equivalent width of zero to those SNe with negative 
measurements. We then selected those SNe with final 
equivalent width errors below 0.6 A. 

For t he photometry, we used d ata from either CfA3 and 
CfA4 (Hickcn ct al.'''2009', "2012^, CSP (jContreras et al.l 
120101 iStritzingcr et al. 2010, 20ll|) and other sources 
(see Table [ij. When multiple sources were available, 
we use either only CSP or only CfA data depending 
on the numbers of photometric points available. We 
corrected for Milky Way extinction using the maps of 
iSchlegel et al.l ()1998D and performed K-corrections and 
S-corrections to the natural rest-frame system of the Las 
Campanas Observatory (LCO) (Stritzingcr ct al. 2011!) 
using photonietry- adiusted spectral template series from 
iHsiao et"aLl (|2007D . In those cases where the natural sys- 
tem of the instrument is not provid ed and t h e phq tom- 
etry is only given in the system of iLandoI^ ()1992D . we 
transform from there to the LCO system. We selected 
only those SNe la with good photometric coverage during 
the Lira law regime (between 35 and 80 days after max- 
imum), which we defined as those SNe having at least 
three simultaneous B and V photometric data points 
with a minimum separation of 25 days between the first 
and last observations. 

For the final sample, we fitted a linear relation in the 
B — V YS time plane between 35 and 80 days after max- 
imum. This led to two values and associated errors for 
every SNe la: the B — V color at 55 days and the Lira 
law B — V decline rate. We retained only those SNe la 
that had an error associated with the Lira law B — V 
decline rate smaller than 0.006 mag day~^. 

For every SN in the sample we al so fitted stretch an d 
colors at maximum using SiFTO (jConlev et al.l 120081) , 
requiring at least one data point between -15 and 7.5 
days after i3-band maximum and another one between 
2.5 and 35 days after _B-band maximum in each band, 
as well as Attt-i.c; values and E(B — V) extinction values 
using SNooPY (jBurns et al.ll20lil ). For most purposes, 
stretch and Ami5 are equivalent quantities describing 
the light curve shape diversity of SNe la, but we use 
both to account for systematic effects in our analysis. 
We also obtained inclinations and host galaxy morpho- 



TABLE 1 

Near by SNe Ia used in this analysis besides data f rom the Center for Astrophysics (CfA3 and Cf A4: IHicken et al.|[2"009I . 
120121) AND THE Carnegie Supernova program ('CSPMContreras ft al.||201(]| : [5tritzinger et al.II201CII) . Not e that we do not 

INCLUDE data FROM THE YALO TELESCOPE BECAUSE OF P OSSIBLE INCONSIS TENCIE S WITH ITS PHOTOMETRIC SYSTEM (IKrISCIUNAS ET AlI 

[2003: .Phillips ft aT1|2006I) . 



Name 



Source 



SN 1993ae 
SN 1994D 
SN 1994ae 
SN 1995D 
SN 1995al 
SN 1995bd 
SN 1996C 
SN 1997E 
SN 1997cw 
SN 1998ab 
SN 1998aq 
SN 1998bp 
SN 1998dm 
SN 1998dh 
SN 1998ef 
SN 1998es 
SN 1999aa 
SN 1999ac 
SN 1999gd 
SN 1999gh 
SN 1999dq 
SN 1999gp 
SN 2000B 
SN 2000CX 
SN 2000dk 
SN 2000fa 
SN 2002aw 
SN 2002bo 
SN 2002CS 
SN 2002dl 
SN 2003cg 
SN 2003du 
SN 2004bg 
SN 2004bk 



Riess et al. (1999) 
Meikle et al. (1996): Altavilla et al. (20041 



Riess et al. (2005* ): Altavilla et al. (2004) 
Riess et al. (1999* ): Sadakane et al. (1996') 

Riess et al. (1999') 

Riess et al. (1999' ): Altavilla et al. (2004) 

Riess et al. (1999') 

Jha et al. (2006) 

Jha et al. (2006') 

Jha et al. (2006') 

Riess et al. (2005') 

Jha et al. (2006) 

Jha et al. (2006) : Ganeshalingam et al. (2010) 

Jha et al. (2006) : Ganeshalingam et al. (2010) 

Jha et al. (2006) : Ganeshalingam et al. (2010) 

Jha et al. (2006) 

Krisciunas et al. (2000) : Altavilla et al. (2004) : Jha et al. (2006) 

Jha et al. (2006) : Phillips et al. (2006) 

Jha et al. (2006) 

Jha et al. (2006 ): Ganeshalingam et al. (2010) 

Jha et al. (2006) 

Jlm^t_^l^ (2006) ; ^aneshalingan^_et_^l^ (2010) ; JCrisciunas_et_aL (20W ) 

Jha et al. (2006) 

2003) : Altavilla et al. (2004) : SoUerman et al. (2004) : Jha et al. 



(2010) 
(2010) 



Candiaetal 

^h^_et_^l^ (2006 ) ; ^aneshalingan^_et_al 

^ha_et_^l^ (2006 ) ; ^aneshalingaii^_et_al 

^aneshalingain^t_^l^ (2010) 

j^risciunas_et_aL (2004 ) ; ^enett^_et_aL (2004) ; ^aneshalingan^_et_al 

^aneshalingain^t_^l^ (2010) 

^aneshalingain^t_^l^ (2010) 

^lias;Ros^^t_aL (2006 ) ; ^aneshalingan^_e^_aL (2010) 

^nu£am^_et_^l^ (2005) ; ^tanishe^^t_ak (2007) ; ^eonar^_et_al^ (2005 ) 

Gan eshalingam et al. (2010) 

Ganeshalingam et al. (2010) 



'2006) 



2010) 



logica l types using the Asiago catalogue (jBarbon et al.l 
Il999f l when available. In order to discard the effects in- 
troduced by possibly very different populations in the SN 
la sample, we have also removed sub-luminous SNe la, 
which we define as those SNe la whose stretch is smaller 
than 0.7 and whose B — V color at maximum is higher 
than 0.4. 

In Figure [T] a) the B — V color evolution for selected 
SNe is shown and in b) for the entire sample used in 
this analysis, separated by Lira law B — V decline rate. 
Figure [T] a) hints that different Lira law B ~ V decline 
rates occur among different SNe la. In the following we 
explore whether the Lira law B — V decline rate is related 
to other observational characteristics in the sample. 



3. RESULTS 

3.1. Lira law B — V decline rates vs equivalent widths 
and colors 

In Figure [2] we show the equivalent widths of blended 
narrow Na I Dl & D2 absorption features vs the Lira 
law B — V decline rate according to the sample defined 
in Section [21 From the Figure it is clear that the distri- 
bution of equivalent widths of those SNe with faster Lira 
law B — V decline rates is different from that of SNe la 
with slower B — V decline rates. If we divide the sample 
by a Lira law B — V decline rate of -0.013 mag day~^, we 
obtain with a Kolmogorov-Smirnov (KS) test that there 
is a probability of 0.0006 that the resulting distributions 
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Fig. 1. — a.) B — V evolution for selected SNe and b) for the entire sample defined in Section (2] dividing the sample between SNe with 
Lira law B — V decline rates faster (red stars) and slower (blue dots) than -0.013 mag day~^. We plot the median evolution and observed 
dispersion for each sample with red and blue lines, respectively. It can be seen that faster Lira law B — V decliners are redder at maximum 
and 55 days afterward. 
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Fig. 2. — a) Equivalent widths of narrow Na I Dl & D2 blended absorption features vs the Lira law B — V decline rate for fast (red stars) 
and slow (blue dots) Lira law B — V decliners and b) the corresponding Kolmogorov Smirnov (KS) test for the cumulative distribution 
of equivalent widths of narrow Na I Dl & D2 blended absorption features for fast (red) and slow (blue) Lira law B — V decline rates for 
the sa mple defined in Section[2] The division between fast and slow Lira law B — V decliners is defined to be at -0.013 mag day~^ (see 
Section las. 



of equivalent widths arise from the same parent popu- 
lation. This division of the sample is based on finding 
the Lira law B — V decline rate that maximizes the dif- 
ferences between fast and slow Lira law B — V decliners 
(explained in more detail in Section 13. 4|) and it implies 
that 33% of the SNe in the sample are defined as fast 
Lira law B — V decliners. 

In order to study the effect of positron absorption, we 
investigate whether the differences between fast and slow 
B — V Lira law decliners are related to their bolometric 
decline rates during the Lira law, whic h we approximate 
by th eir F-band decline rates (see e.g. iCappellaro et al.l 
Il997f ) . We found that fast and slow Lira law B — V declin- 
ers do not have significantly different distribution of V- 
band decline rates, with a KS test probability of 0.86 that 
their distributions of y~band decline rates arise from the 
same parent population. However, we find a significant 
difference between their distributions of _B-band decline 
rates, with a KS test probability of only 0.0001 that their 
distributions of i3-band decline rates arise from the same 
parent population. Thus, the Lira law B — V decline rate 
diversity is driven primarily by different B-band decline 
rates, in the direction that fast Lira law B — V decliners 
are slow Lira law B declinerfl The fact that the Lira law 
B — V decline rate is connected to the presence of nar- 
row absorption features suggests that this effect is either 
caused by younger progenitor systems having slower Lira 
law B decline rates, which would be associated with star 
forming regions and larger gas/dust column densities; or 
to the presence of circumstellar material that causes the 
B band evolution to become slower at late times. We 
show the B and V band decline rates vs the B — V de- 
cline rates in Figure [3] 

If faster Lira law B — V decliners are associated with 
more material in the line of sight, we should also expect 
a difference in the distributions of colors around maxi- 
mum light of fast and slow Lira decliners. In fact, the 
distributions oi B — V colors at maximum obtained with 
SiFTO and E{B-V) extinctions obtained from SNooPY 
are significantly different if we divide the sample between 
fast and slow Lira law B — V decliners, with p-values of 
0.0004 and 0.0001, respectively. This is shown in Fig- 
ure 01 Fast and slow Lira law B — V decliners also differ 
significantly in their colors 55 days after maximum {p- 
value of 0.003). It should be noted that fast and slow 
Lira law B — V decliners tend to have more similar col- 
ors at late times during the Lira law, as redder SNe in 
B — V tend to decrease faster \n B — V . These results 
are summarized in Tabled 

3.2. Environmental and SN intrinsic effects 

After finding that differences in the measured Lira law 
B — V decline rates are related to the presence of mate- 
rial in the line of sight, we tested whether this could be 
explained by faster Lira law B — V decliners occurring in 
younger stellar environments, i.e. whether the material 
detected could be ISM. In order to do this we divide our 
sample by quantities that should favor higher ISM col- 
umn densities, such as host galaxy morphological type or 
host galaxy inclinations (removing E/SOs from the sam- 

^ Note that this is different to the _B-band decHne rate directly 
after maximum, which is connected to different stretch or Amis 
values 



pie, where it is difficult to define an inclination). If there 
was a strong preference for morphological type among 
different Lira law B — V decline rates it could mean that 
younger systems decay faster during the Lira law, and 
that there is no need to invoke CSM as an explanation. 
If there was a strong preference for galaxy inclination, 
it would suggest that the Lira law B — V decline rate is 
affected by ISM in the line of sight. 

Stretch and Ami5 are also thought to be as- 
sociated w ith the age o f the p rogenitor systems 
(iHainuv et al. 1995. 1996; Ri esset"aLll l999: Hamuy et aD 
l2000t iSullivan et al.l [20061 ). which we also confirmed in 
our sample, withp-values of less than 0.00001 for the null 
hypothesis that the distributions of stretch and Amis 
dividing the sample between earlier and later type hosts 
arise from the same parent population. Thus, stretch and 
Amis pla-y a similar role as the host galaxy morphologi- 
cal type in the sense that younger systems (larger stretch 
or smaller A77115) occur in regions where the ISM column 
density should be larger. Therefore, if stretch or Amis 
were connected to the evolution of B — V at late times it 
could suggest that environment could drive the different 
Lira law B — V decline rates. On the other hand, stretch 
and Am, IS are known to account for most of the diversity 
of SN la light curves because they are likely related to 
the enclosed mass of ^^Ni or iron group elements (IGE) 
in the ejecta, which is the power source of SN la light 
curves. Considering that at late times the physics of the 
decay of ^^Co and the absorption of positrons start to 
become more relevant, it is not unreasonable to think 
that both stretch and Am,is could be related to the Lira 
law B — V decline rate. 

Thus, we test whether the distributions of host mor- 
phological type, host inclination, stretch and Attiis of 
fast and slow Lira law B — V decliners arise from the 
same parent population. In none of the tests we obtained 
strongly significant differences (see Table [3]). Therefore, 
assuming that the Lira law B — V decline rate is not 
strongly affected by the amount of ISM in the line of 
sight, we interpret the results shown in Table [2] as evi- 
dence for CSM around those SNe la having the fastest 
Lira law B — V decline rates (but see Section IX6|) . 

The presence of CSM could easily change the slope 
of the Lira law B — V decline rate, either by mul- 
tiple scattering processes ( i.e. CSM light echoes, see 
lAmanullah fc Goobad 120111 ) or by the progressive de- 
struction of nearby dust if the material is sufficiently 
close to the expanding ejecta, making the B — V colors in- 
creasingly bluer at late times and, as a consequence, the 
Lir a law B — V d e cline r ates faster. This was also noted 
by iFolatelli et al.l ()2010D in systems with color excesses 
that were moderate to large. 

If the B — V colors at maximum are mainly driven by 
CSM reddening and if either light echoes or late time 
CSM dust destruction occur, we would expect that the 
reddening observed during maximum light should de- 
crease at later times, making the reddening to be first 
dominated by CSM and later by ISM. In fact, we cannot 
rule out that the the distributions of colors at maximum 
in the earlier and later type hosts - or lower and higher 
inclination hosts - arise from the same parent population 
(p-values of 0.16 and 0.39, r espectively), as also f ound 
in previous studies (see e.g. [Johansson et al.ll2012l c.f. 
ISullivan eFal..,2010, ). This would be expected if the col- 



6 Forster, Gonzalez-Gaitan, Folatelli & Morrell 

TABLE 2 

KS TEST FOR THE DISTRIBUTION OF OBSERVED PROPERTIES DIVIDING THE SAMPLE BY LiRA LAW B — V DECLINE RATES. THE NULL 

HYPOTHESIS IS THAT FAST AND SLOW LiRA LAW B — V DECLINING SNe Ia HAVE THE SAME DISTRIBUTION OF BLENDED NA I D1 & D2 

NARROW ABSORPTION FEATURES, COLORS AT MAXIMUM (SiFTO), EXTINCTION VALUES (SNOOPY) OR COLORS 55 DAYS AFTER MAXIMUM. 



KS test EW(Na) B - V a.t max. E{B - V) B - V at 55 days 



p-valuc 



0.0006 



0.0004 



0.0001 



0.003 



TABLE 3 

KS TEST FOR THE DISTRIBUTION OF OBSERVABLES DIVIDING THE SAMPLE BY LiRA LAW B — V DECLINE RATES. THE NULL HYPOTHESIS IS 

THAT FAST AND SLOW DECLINING LiRA LAW SNe IA HAVE THE SAME DISTRIBUTION OF HOST MORPHOLOGICAL TYPES, HOST INCLINATIONS, 

STRETCH OR Amis. THE NULL HYPOTHESIS CANNOT BE DISCARDED IN ANY OF THESE TESTS, SUGGESTING THAT THE EQUIVALENT WIDTH 

AND COLOR DIFFERENCES SEEN IN FIGURE [2] ARE NOT DRIVEN BY ENVIRONMENTAL FACTORS. 



KS test 


Host morph. type 


Host. incl. 


stretch 


Amis 


p-value 


0.10 


0.19 


0.18 


0.07 



ors of SN la at maximum were primarily driven by ISM. 
However, we do rule out that the distributions of colors 
55 days after maximum in earlier and later type hosts 
arise from the same parent population (p-value of 0.003). 
This suggests that: 1) B — V colors at maximum and 55 
days after maximum are driven by different physical ef- 
fects, and 2) B — V colors 55 days after maximum are 
more related to the properties of the environment where 
SNe occur. This is consistent with the CSM interpreta- 
tion presented above when the Lira law decline rates are 
considered. These results are summarized in Table U) 

We should emphasize that we are not suggesting that 
ISM absorption is irrelevant at maximum light, but only 
that it does not drive the results shown in Table [2] (al- 
though see Section lS^ . In fact, we detect a highly signif- 
icant increase on the equivalent widths of Na I Dl & D2 
due to ISM related properties by performing KS tests on 
the distribution of EWs of Na I Dl & D2 dividing a big- 
ger sample by host galaxy type, host galaxy inclination, 
stretch or Amis instead of by Lira law B — V decline 
rate, with p~values smaller than 0.00001 in all cases. 

3.3. V — i colors and dust reddening laws (Rv) 

A prediction from IGoobari (|2008D is that those SNe 
la with more CSM should have a steeper reddening law 
towards the blue. Thus, we expect that fast Lira law 
B — V decliners should have a lower Rv at maximum 
light if CSM was present in these systems, which can be 
measured using V — i colors. In Figure [5] we show the 
V — i vs B — V colors at maximum and 55 days after max- 
imum with best-fitting Ry valu es ass uming the r edden - 
ing law of lCardenietall (IT989I) and lO'Donnelll (Il99l . 
To estimate the Rv best-fitting errors we have done a 
Montecarlo analysis based on the measured B, V and 
i magnitudes and their errors. We created 1000 real- 
izations of these magnitudes at maximum and 55 days 
after maximum to derive B — V and V — i colors, which 
by definition are correlated, and then derive 1000 best- 
fitting i?v values, whose standard deviation was used as 
the error. 

We find that fast Lira law decliners have a steeper 
reddening law (lower Rv) than slow Lira law declin- 
ers, a result that becomes more pronounced at max- 
imum light, even after removing the highly reddened 
SN 2003cg and SN 2006X. In fact, the i?v values that 
best match the data for fast and slow Lira law decliners 



are i?v = 1-90 ± 0.09 and i?v = 3.23 ± 0.29, respec- 
tively, which is consistent with wha t has been found for 
highly reddened and normal SNe la (jFolatelli et al.lf2010l : 
Mandel et al. 2011" ). At 55 days after maximum we find 
that the reddening laws of fast and slow Lira law declin- 
ers are significantly different, but more similar to each 
other than at maximum light, with Ry = 2.61 ± 0.05 
and Ry = 3.63 ± 0.12 for fast and slow Lira law B — V 
decliners, respectively. These results are consistent with 
a contribution from light echoes and/or CSM dust de- 
struction at late times, which would make the B — V 
colors bluer. However, we cannot rule out an intrinsic 
relation between V ~ i and B — V colors at maximum 
light and at late times that could mimic some of our re- 
sults, although this relation would need to be different 
for fast and slow Lira law B — V decliners. 

We have also studied the V — i colors at maximum and 
55 days after maximum of the SNe in the sample, in a 
similar fashion to B — V colors. The probabilities that 
the distributions oi V — i colors at maximum light or 55 
days after maximum of fast and slow Lira law B — V 
decliners arise from the same parent population is 0.13 
and 0.07, respectively. The lower significance obtained 
with these KS tests may indicate that the effects we are 
measuring are stronger at bluer wavelengths. 

3.4. Robustness and systematic errors 

In the previous analysis we set a seemingly arbitrary 
division of the sample at a Lira law B — V decline rate 
of —0.013 mag day~^. In order to define this limiting 
decline rate we performed a series of KS tests for the dis- 
tributions of equivalent widths of blended Na I Dl & D2 
narrow absorption and B — V colors. The optimal divi- 
sion was defined by varying the dividing Lira law decline 
rate and finding a small combined p-values of the KS- 
tests for the distributions of equivalent widths and B — V 
colors, which corresponded to a dividing Lira law B — V 
decline rate of about -0.013 mag day~^. This is shown 
in Figure [6l where we also show the p-values obtained 
for the KS test associated with some of the environmen- 
tal variables discussed before. This Figure shows that: 
1) the equivalent width of blended Na I Dl & D2 nar- 
row absorption features and the colors at maximum and 
55 days after maximum have a minimum p-value at a 
similar dividing Lira law decline rate and 2) the envi- 
ronmental variables considered in this analysis are never 



TABLE 4 

KS TEST FOR THE DISTRIBUTION OF B — V COLORS AT MAXIMUM (B — V max) AND 55 DAYS AFTER MAXIMUM [B — V 55), DIVIDING THE 

SAMPLE BY Lira law B — V decline rates (columns 1 AND 2) OR BY THE MEDIAN MORPHOLOGICAL TYPE (COLUMNS 3 AND 4). THE 

NULL HYPOTHESIS IS THAT FAST AND SLOW LiRA LAW B — V DECLINING SNe IA HAVE THE SAME DISTRIBUTION OF COLORS AT MAXIMUM 

LIGHT AND 55 DAYS AFTER MAXIMUM LIGHT, OR THAT EARLIER AND LATER TYPE HOSTS HAVE THE SAME DISTRIBUTION OF COLORS AT 

MAXIMUM AND 55 DAYS AFTER MAXIMUM LIGHT. THESE RESULTS SUGGEST THAT THE REDDENING EFFECT OF CSM IS MORE APPARENT AT 

EARLY TIMES AND THAT OF ISM, AT LATE TIMES. 
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strongly associated with different Lira law decline rates, 
independently of the division of the sample used in the 
KS tests. Comparing the range of Lira law decline rates 
considered with those found in Figure [U we see that the 
high significance is not very sensitive to the dividing de- 
cline rate, making our conclusions more robust. 

In order to check whether the results shown in the pre- 
vious sections could be an artifact due to different mea- 
surement errors, we performed a KS~test on the distribu- 
tion of the errors associated with the Lira law B — V de- 
cline rate among the fast and slow Lira law B — V declin- 
ers, obtaining no significant differences between them (p- 
value = 0.34). In addition, in order to test whether our 
conclusions could be affected by constraining the sample 
to SNe with equivalent width errors lower than 0.6 A and 
to SNe with Lira law B — V decline rate errors lower than 
0.006 mag day~^, we relaxed both constraints doubling 
the required errors and obtained consistent results in all 
the relevant tests. Furthermore, to test whether fast and 
slow Lira law decliners are distributed differently with 
distance, which could affect some of our results through 
K-corrections, we performed a KS-test on the distribu- 
tion of the SN distances separating the sample between 
fast and slow decliners, obtaining no significant differ- 
ences among the two groups (p-value — 0.49). To test 
this possible bias differently, we restricted our sample to 
SNe with distance modulus smaller than 35 mag and re- 
peated all the KS-test from previous sections, obtaining 
consistent results again. Restricting the sample to only 
CfA3, CfA4 and CSP photometry we also recover our 
main results, including different reddening laws between 
fast and slow Lira law B — V decliners. 

We have also repeated our KS tests above after remov- 
ing those SNe la that are known to be extreme cases of 
either high column densities, variable Na lines or peculiar 
ejecta properties: SN 2003cg, SN 2006X, SN 20071e, SN 
2000cx and SN 2005hk, and our main conclusions remain 
the same. This indicates that the correlations found ap- 
ply to a high fraction of SNe la, and not only to those 
SNe la which have strong CSM-like features like variable 
Na absorption. 

Finally, in order to test whether the time interval used 
for measuring the Lira law B — V decline rate (35 to 
80 days after maximum) affects our conclusions, we have 
repeated the relevant tests using data from different time 
intervals after maximum. Later time photometry could 
have larger measurement errors and could therefore affect 
our best fitting decline rates. Also, the choice of the start 
of the Lira law regime is important since it could be too 
close to the observed peak in color evolution. Thus, we 
measured decline rates with the following time intervals 
after maximum: 30 to 90 days, 40 to 80 days, 30 to 



70 days and 50 to 90 days. In all but the last interval, 
between 50 to 90 days, we recover our main results. This 
could be due to the effect of larger errors or fewer data 
points at later times, or alternatively, to a systematic 
flattening of the decline rate at later times in the fast 
Lira law decliners. 

3.5. Nebular velocities 

In FG12, a correlation between the presence of narrow 
absorption features and positive nebular velocity shifts 
(uncb, introduced bv lMaeda et al.l 20101 was found, which 
was interpreted as evidence for asymmetries in the CSM 
post explosion. In order to test whether the correlations 
found by FG12 are consistent with the results presented 
in this work, in Figure [7] we show the nebular velocity 
shift of nearby SNe la against their Lira law B — V de- 
cline rate. For this we have taken a we ighted average 
of the nebular velociti es from iBlondin et al. (2012) and 
[Silverman et al.l(j2013[ ). We could not reject the hypothe- 
sis that the distributions of nebular velocities of fast and 
slow Lira law B — V decliners differ (p-value of 0.09). 
Nonetheless, the number of SNe in the test is small and 
by changing the definition interval for the Lira law B — V 
decline rate from 35 to 80 days to 30 to 70 days we ob- 
tain a p-value of 0.003. This may suggest that SNe la 
with the fastest Lira law B — V decline rates tend to 
have positive nebular velocity shifts. These are also the 
SNe found to have stronger Na I Dl & D2 absorption in 
FG12. 

We note that in Figure [7] the fast declining SN with 
Wncb < is SN 1994D, which is the only supernova in 
the sample whose Lira law B — V decline rate differs sig- 
nificantly from the subtraction of its Lira law B and V 
decline rates, used in Figure [3l If the latter is used, its 
decline rate becomes much lower, -0.011 ± 0.002 mag 
day^^ instead of -0.015 ± 0.001, a difference which is 
explained by non-simultaneous B and V photometry for 
some days, which cannot be used to measure B — V de- 
cline rates directly. If this different value is used, we 
obtain a p-value of 0.02 for the probability that the dis- 
tributions of nebular velocities of fast and slow Lira law 
B — V decliners arise from the s ame population. 

Given the correlation found bv lMaeda et al.l (J2010D be- 
tween nebular velocities and broad Si II A6355A absorp- 
tion velocity gradients, we investigated whether the Lira 
law B — V decline rates are related to the Si II absorp- 
tion velocity at maximum. We measured the Si II ve- 
locities at different epochs and obtained values at maxi- 
mum by interpolation, then combin ed these values with 
measurements from other sources (iBenetti et al.l 120051: 
Foley et al.l I20TTI: iSilverman et~all 120121: iFolatelli et all 
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Fig. 3. — a) Lira law B-band decline rates vs Lira law B — V decline rates for fast (red stars) and slow (blue dots) Lira law B ~ V 
decliners and b) Lira law V-band decline rates vs Lira law B — V decline rates for fast (red stars) and slow (blue dots) Lira B ~ V law 
decliners. We found that the diversity of iJ — V decline rates is mostly due to different B-band decline rates during the Lira law, and not 
to V— band decline rates variations, suggesting that the origin of the diversity has a larger effect at shorter wavelengths. 
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Fig. 4. — a) KS tests for the distributions of B — V colors at maximum obtained from SiFTO dividing the sample by the Lira law B — V 
decline rate and b) KS tests for the distributions of E(B — V) extinction values obtained from SNooPy dividing the sample by the Lira 
law B — V decline rate. Fast Lira law B — V decliners are r edder at maximum, have stronger Na I Dl & D2 absorption, but have no strong 
preference for any particular environment (see Section l3.2l l. which is consistent with them having more CSM. 
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Fig. 5. — a) V — i and B — V color at maximum light and b) 55 days after maximum light, dividing the sample between fast (red stars) 
and slow (blue dots) Lira law B — V decliners. We fit a straight line to the Y — i vs B — V colors at maximum light and 55 days after 
maximum for each sample. At maximum light and 55 days after maximum the extinction laws are significantly different, with significantly 
lower Ky for fast Lira law B — V decliners at all times, although becoming more similar to slow Lira law B — V decliners 55 days after 
maximum light. We show for reference extinction vectors with Ay = 1 for Ry = 3.2 (continuous) and Ry = 1.7 (dashed). If we remove 
the highest extinction SNe from the fast declining sample we obtain an even lower Ry at maximum. 



rors to obtain averages. We tested whether Si II ab- 
sorption velocities at maximum are related to Lira law 
B — V dechne ra tes, c olors at ma ximum (as found by 
iWang et alJ 120091 and iFolev et all [2011) or equivalent 
widths of blended Na narrow absorption features, which 
would be consist ent with the excess blueshift found by 
iFolev et"aLll2012l We found that there is a probability 
of 0.07 that the distributions of Lira law B — V decline 
rates of SNe la with Si II velocities faster or slower than 
11,800 [km s~^] arise from the same parent population, 
a probability of 0.002 that the distributions of colors at 
maximum of fast and slow Si II velocity SNe la arise 
from the same parent population, a probability of 0.81 
that the distributions of equivalent widths of narrow Na 
absorption of fast and slow Si II velocity SNe la arise 
from the same parent population and a probability of 
0.07 that the distributions of morphological types of fast 
and slow Si II velocity SNe la arise from the same par- 
ent population. This confirms the color differences found 
bv lWang et al.l ()2009[ ) and suggests a further connection 
with Lira law B — V decline rates. The non-detection 
of significant differences in the distribution of equivalent 
widths of Na absorption or morphological type with Si II 
velocity may indicate that some of these color differences 
are intrinsic, or that our sar nple size is not larg e enough, 
given the recent finding by (|Wang et al.ll2013f ) of differ- 
ent radial distributions for fast and slow Si II velocity 
SNe la. 

3.6. ISM as the origin of Lira law decline rate 
variations? 

Although the previous tests suggest that environmen- 
tal factors do not contribute to the Lira law B — V decline 
rate, it is intriguing that in Figure[6]the region where the 
p-values associated with tests for a relation between Lira 
law decline rates and the host galaxy morphological type 
or inclination are lowest coincides with the region where 
p-values associated with equivalent widths or colors are 
also the lowest. This led to us to consider whether the 
ISM could really affect the Lira law B — V decline rate 
in some way. 

Since the inclination of a galaxy should not be related 
to the presence or absence of CSM in SNe la, and given 
that the ISM will be generally too distant from the SN 
explosion to produce significant light echoes or to be 
destroyed by the expanding ejecta between 35 and 80 
days after maximum, a different mechanism to produce 
changes in the Lira law B — V decline rate would be 
needed. Note that we ar e not referring to ISM light 
echoes (|Wang et al.l I2008D , which can affect the light 
curves at even later times than those considered in this 
work. In fact, the transversal expansion of the ejecta in 
an inhomogeneous distribution of projected column den- 
sities of the host galaxy ISM could affect the Lira law 
B — V decline rate by changing the average column den- 
sity as the photospheric radius increases. If this is the 
case, what needs to be shown is whether systems that 
have stronger Na I Dl & D2 absorption and redder col- 
ors at maximum are also expected to have faster Lira law 
B — V decline rates due to this effect^ 

For this, we follow the approach of iPatat et al.l (|20100 
to study the effect of the transversal expansion of the 
ejecta in a non-homogeneous ISM. Assuming their ap- 
proximation for the evolution of the SN photospheric ra- 



dius with time, and using random realizations of the ISM 
column density assuming a fractal structure with power 
law of -2.75 on spatial scales down to the size of the pho- 
tosphere, we compute the average column density within 
the photospheric radius as the SN expands. In Figure [51 
left, we show one realization of the dust column density 
with different photospheric radii shown as the SN ex- 
pands. In this example, because the SN started behind a 
relatively over-dense region, the average column density 
started at a relatively large value and tended to decrease 
with time due to a dilution effect as the photospheric ra- 
dius increases. In Figure [3 right, we show the expected 
evolution of the average column density for one thou- 
sand ISM realizations, showing different percentiles for 
reference. This figure shows that the biggest dispersion 
occurs at early times, when chance alignments with un- 
der or over-dense regions can dominate, but the dilution 
effect occurring at later times will make the dispersion 
of the average column densities decrease. 

Thus, Figure[8|shows that the inhomogeneous nature of 
the ISM will naturally produce SNe la that have stronger 
equivalent widths of narrow blended Na I Dl & D2 ab- 
sorption and are redder at maximum, and that have fast- 
evolving average column densities, and therefore faster 
Lira law B — V decline rates. The higher the inclina- 
tion, the stronger the effect as the typical column density 
will tend to increase. Moreover, at later times the colors 
will be dominated by the overall column density, which 
could explain why SNe la in later morphological types 
are significantly redder at later times. However, this ef- 
fect would not explain the differences in reddening laws 
observed between fast and slow Lira law B — V decliners. 

Nonetheless, the column density is expected to be pro- 
portional to the product of the color excess and Ry, and 
a detailed analysis of Figure [8| implies that at maximum 
light the B — V colors of the SNe in the percentiles 50 
and 95 should differ by about 0.25 ma^, or by about 
0.12 mag for the percentiles 30 and 70. At 55 days af- 
ter maximum, the B — V colors of the SNe in the per- 
centiles 50 and 95 should differ by about 0.15 mag, and 
by about 0.08 for the percentiles 30 and 70. These mag- 
nitude differences are of the same order as those seen in 
Figure [TJ but are not consistent with the lack of corre- 
lation between the Lira law F-band and B — V decline 
rates observed in Figure [S}). This suggests that the in- 
homogeneous nature of the ISM can not explain all our 
findings, and that a more detailed analysis will be needed 
to exclude this possibility. 

4. DISCUSSION AND CONCLUSIONS 

In this work we have shown that there are significant 
differences in the B — V evolution of SNe la during the 
Lira law regime, between 35 and 80 days after maximum. 
These differences are driven by the evolution of SNe la in 
the _B-band: fast Lira law B — V decliners are slow Lira 
law B decliners. With the sample defined in Section [2] we 
have shown that fast Lira law B — V decliners, defined 
as having Lira law B — V decline rates smaller than - 
0.013 mag day~^, have stronger narrow absorption lines 
of blended Na I Dl & D2, are redder and have a a lower 
Ry extinction law at maximum light and 55 days after 



3 -2.51og(^)/iJv = -0.25 assuming Ry 
units and nso = 0.9 units. 
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Fig. 6. — KS test p— values for the null hypothesis that the distribution of the variables indicated by the labels arise from the same parent 
population when dividing the sample by the Lira law B — V decline rate shown in the a;-axis. This plot shows that the distribution of 
colors at maximum and equivalent widths of blended Na I Dl &; D2 narrow absorption features have a robust correlation with the Lira law 
B — V decline rate, and that the division of the sample between fast and slow Lira law decliners by a Lira law B — V decline rate of -0.013 
[mag day^^] is appropriate. 
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Fig. 7. — Nebular velocity shift of SNe la against their Lira law B — V decline rate. Faster Lira law B — V decliners appear to have 
positive nebular velocity shifts. The fast declining SN with Dncb < is SN 1994D (see discussion in the text). 
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maximuni light. The differences in colors were stronger 
at maximum light. 

Additionally, we did not find a highly significant dif- 
ference between the distributions of host galaxy morpho- 
logical type, host galaxy inclination, stretch or Attiis di- 
viding the sample by Lira law B — V decline rate. This 
suggests that the Lira law B — V decline rates are not 
strongly dependent on the age of the SN la progenitor 
systems, are not strongly affected by ISM, and are not 
mainly driven by the the mass of nickel in the ejecta, the 
main parameter of SN la light curves. 

For the sample of 89 SNe la that had good photometric 
coverage up to late times, we found that their colors at 
maximum were consistent between earlier and later type 
galaxies. The colors 55 days after maximum were sig- 
nificantly different between earlier and later type galax- 
ies, suggesting that the extinction at later times is more 
clearly connected to ISM and that late time colors may 
be a better proxy for environment. If CSM is dominant 
at early times, the cancelling effect on the color excess 
due to CSM light echoes or to dust destruction as the 
ejecta sweeps over the CSM could explain this result. 

A SN explosion surrounded by CSM should sublimate 
dust and ionize Na I up to some maximum radii, there- 
fore the CSM should have a characteristic minimum ra- 
dius, rcsM, which is thought to be between 0.01 an d 0.1 
pc for Na I (jPatat et al.ll2007t iSimon et al.ll2009D . and 
which is expected to be of the same o rder of magnitude 
for dust (AmanuU ah fc Goobadl2011f l. Photons leaving 
the ejecta at maximum light could be seen as light echoes 
55 days after maximum if they are reflected after travel- 
ing for about 55 days, or a typical distance of about 0.05 
pc, which is of the order of this minimum CSM radius. 
However, light echoes are not the only mechanism that 
can produce this effect: depending on the outer ejecta 
speed, which we assume to be one tenth of the speed of 
light, the ejecta could have traveled up to about 0.001 pc 
at maximum light; therefore it should not have time to 
interact with the CSM and sublimate CSM dust. How- 
ever, 55 days after maximum light the ejecta could be 
approaching the inner edge of the CSM, since it will have 
traveled up to about 0.006 pc, and it could start reducing 
the dust column density. This could explain why the dif- 
ferences in colors and extinction laws found at maximum 
light between fast and slow Lira law decliners were not 
as strong 55 days after maximum, but also why environ- 
mental effects explain some of the differences in colors 55 
days after maximum. 

It is interesting to note that smaller dust sublimation 
or Na I Dl & D2 ionization CSM radii should produce 
more extinction and stronger narrow absorption features, 
which we find in fast Lira law B — V declining SNe la. 
For the same ejection velocity and mass loss rate be- 
fore a supernova explosion, the dust column density and 
extinction should scale as T'^gj^.j- Smaller CSM dust in- 
ner radii should also produce stronger CSM light echoes, 
since their intensity should scale as r, 



CSM ' 



which would 
be reflected in faster B — V evolution at late times be- 
cause of the slower expected Lira law i3- band decline 
rates. In fact, lAmanullah fc Goobaii (|2011f) showed that 
smaller CSM radii should affect more significantly the 
B — V evolution at later times, producing faster Lira law 
B — V decline rates. 



The resulting fast declining sample is composed of 
29 SNe la (33%) and the slow declining sample, by 60 
SNe la (67%). Interestingly, the fraction of fast Lira 
decliners is similar to the fraction of 35% high veloc- 
ity SNe la (HV SNe la, with ^;(S'i//)max > 11,800 



[km 



I k m s I) amo n g tne group ot normal blNe ia louna 
lay iWang et all (|2009| ). or the excess of about 25% of 
SNe la with blueshifte d absorption features detected by 
iSternberg et al.l (j2011h . which is also consistent with the 
excess of blueshifted systems among HV SNe la found by 
[Fblcy et al . (2012). This suggests that blueshifted nar- 
row absorption features, high velocity Si II features at 
maximum, high equivalent widths of narrow absorption 
features, redder colors at maximum, fast Lira law decline 
rates and lower Ry reddening laws may all be associated 
with the same population, which in our interpretation 
is the one with significant CSM. Unfortunately, there 
is no overlap between our sample and the recent sam- 
ple of CSM-interacting SNe la found bv [Silverman et ahl 
(l2013fl. The overlap b etween our sample and that from 
ISternberg et al.l (120111 ) is very small (only 4 SNe) and 
we cannot test whether their blueshifted sample corre- 
sponds to our fast declining sample. However, the two 
overlapping SNe la in their blueshifted sample are also 
fast Lira law B ~ V decliners in our sample, SN 2006X 
and SN 20071e, and they also happen to h ave variable 
Na I pi fc D2 narrow absorption features (jPatat et al.l 
I2007t ISimon et al. 2009), which is expected in the CSM 
interpretation. 

In the context of FG12, we have found a weak correla- 
tion between nebular velocity shifts (fneb) and Lira law 
B — V decline rates, which is consistent with the inter- 
pretation of FG12, but which needs to be confirmed with 
a bigger sample0. This would mean that SNe la with 
"yncb > have stronger equivalent widths of Na I Dl & 
D2, are redder at maximum light, have faster Lira law 
B — V decline rates and have lower Ry reddening laws at 
maximum light because they have more CSM in their line 
of sight. Assuming that about 33% of the SNe la have 
significant CSM in the line sight and that Wnob is a good 
proxy for explosion viewing angle, this would mean that 
there is a comparable population of SNe la with negative 
nebular velocities that would also have significant CSM, 
but which we do not see in narrow absorption features 
or fast Lira law B — V decline rates because their CSM 
is on the opposite side of their ejecta. Here it should be 
not ed that dust scattering is pr imarily forward scatter- 
ing (jHenvev fc Greensteinlll94ll ) , which means that light 



echoes would be best seen when the CSM is between the 
SNe and the observer (wneb > in FG12). Although the 
fractions are uncertain at the moment, this leaves room 
for additional progenitor scenarios that do not fit in this 
CSM picture. 

Since changes in the Lira law decline rate could be eas- 
ily produced by CSM (jAmanullah fc GoobaJl201lL see 
their Figure 8), and given that we did not find any sig- 
nificant environmental dependency of the Lira law B — V 
decline rates, we propose that the differences found be- 
tween fast and slow Lira law B — V decliners are due to 
CSM, but we also investigate the possibility that a sim- 

* We have confimied the correlation between narrow absorp- 
tion Na lines and nebu lar velocities using the public spectra from 
ISilvernian et all l|2013l ) and the method described in FG12. 
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Fig. 8. — Left: example dust realization from the observer's perspective following the procedure of lPatat et al.l l|2010l ) compared to the 
typical size of the SN la photosphere up to 100 days after explosion. Right: Evolution of the average column density within the photospheric 
disc of SNe la in 1000 ISM realizations. Most of the dust column density variations are found at earlier times. This means that the most 
reddened objects at earlier times will be less reddened at later times, implying a faster Lira law B — V decline rate. The dashed— blue line 
corresponds to the example ISM realization shown, which is very atypical among the simulated ISM realizations. 
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ilar result could be produced by ISM. We showed that 
the transversal expansion of the SN photosphere within 
an inhomogeneous distribution of column densities could 
naturally lead to faster Lira law decliners having stronger 
absorption at maximum as discussed in Section [3.61 but 
further modeling is required to exclude this possibility. 

This work is an addition to the growing evidence 
for CSM in SNe la, but we believe that only the 
use of high reso lution spectra (jSternberg et al.l 120111 : 
iFolev et al.l 120121 ) will help better constrain the na- 
ture of the absorbing material in SNe la. The lack 
of evidence for shocked winds, shocked companion 
stars, companion stars in pre-supernova images o r com 
panion stars in SN la remnants (e.g. | Li et al.l 2011 



Nugent et al."2011':'Bloom et al. i20TalBrown et al.ll2012 : 
Schacfcr & Pagnotta 2012; Edw ards et al.ll2012[ ) inav in- 
dicate that the progenitors of SNe la arise from more 
than one progenitor scenario or from more complicated 
scenarios that we have not yet imagined. 

Finally, although it is still possible that the physics of 
positron absorption could be responsible for the varia- 
tions of the Lira law detected, this effect should depend 
on the age of the progenitors in order to produce stronger 
narrow absorption features due to larger ISM column 
densities. The fact that we do not find a strong depen- 
dency between Lira law decline rates and environment, 
stretch or Amis suggests that this is not the case. Ad- 



ditional studies of the relation between very late B — V 
decline rates (e.g. beyond 200 days) and the Lira law 
B — V decline rate could be used to provide a clean sep- 
aration between the effects of CSM and the physics of 
positron absorption (see e.g. ILair et al.|[2006|) . 
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